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PATENT 

Attorney Docket No.: 201 74-0024 lOUS 
Client Reference No.: CIT-3366 

Velocity Independent Analyte Characterization 

STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH 

The U.S. Government has a paid-up license in this invention and the right 
5 in limited circumstances to require the patent owner to license others on reasonable terms 
as provided for by the terms of Grant No. HG-01 642-02, awarded by the National 
Institute of Health. 

CROSS-REFERENCE TO RELATED APPLICATIONS 
This application claims the benefit of U.S. Provisional Apphcation No. 
10 60/237,937, filed October 3, 2000, which is incorporated herein by reference in its 
entirety. 

FIELD OF THE INVENTION 
The present invention relates to an apparatus and a method for determining 
a characteristic parameter of an analyte in a fluid medium. In one particular aspect, the 
15 present invention relates to a velocity independent flow cytometry. 

BACKGROUND OF THE INVENTION 
Recently, there has been a growing interest in microfluidic flow cytometry. 
Several chip based systems have been demonstrated for cytometry and sorting of cells 
and molecules. See, for example, Unger et al., Science, 2000, 288, 1 13-116; Fu et al, 

20 Nature Biotechnol, 1999, 17, 1109-1111; Quake et al., Science, 2000, 290, 1536-1540; 
Schrum et al, Anal. Chem., 1999, 77, 4173-4177; Knight et al., Phys Rev. Lett., 1998, 80, 
3863-3866; Chou et al, Proc. Natl. Acad. Sci. USA, 1999, 96, 1 1-13; Chou et al, 
Electrophoresis, 2000, 21, 81-90; and Kameoka, et al.. Sensors and Actuators B, 2001, 
77, 632-637. There are generally two ways to pump fluid in these devices: pressure 

25 driven flow or by electroosmotic forces. Pressure driven flow results in Poiseulle flow, 
which has a parabohc velocity distribution in the channel. This complicates measurement 
of analytes since each analyte, e.g., cell or molecule, passes through the interrogation 
region with a different velocity. One can mitigate the effects of the Poiseulle flow by 
using a sheath fluid for hydrodynamic focusing, but this introduces other issues by 

30 diluting the sample and complicating downstream analysis. Although electroosmotic 



flow is typically more uniform and plug-like than pressure driven flow, it too results in 
variability in flow velocity. See, for example, Schrum et al.. Anal. Chem., 1999, 71, 
A\l?>-A\ll. In addition, most, if not all, electroosmostic flow requires careful balancing 
of the ions in the solution and attention to prevent ion depletion. Furthermore, in some 
5 cases, it has also been shown that eukaryotic cells are difficult to manipulate 

electroosmotically. See, for example, Li et al, Anal. Chem., 1997, 69, 1564-1568. 

One possible method for measuring a velocity independent characteristic 
parameter, e.g., fluorescence, of an analyte is to use a uniform detection zone, e.g., 
excitation region, large enough to illuminate the entire particle or molecule of interest. In 

1 0 this case, the height of the detected fluorescent peak will be substantially proportional to 
the fluorescence intensity of the particle. See, for example, Chou et al., Proc. Natl. Acad. 
Sci. USA, 1999, P5, 11-13. While this method is substantially not affected by the 
distribution of velocities, using only one point from the entire peak, namely its maximum, 
exploits only a small part of the information that is embedded in the peak. Moreover, the 

1 5 accuracy of this method is susceptible to noises. On the other hand, measuring the area 
underneath the entire peak, which is proportional to total fluorescence intensity integrated 
over the excitation duration, is velocity dependent measiirement. This is due to the fact 
that faster particles have narrower peaks than slower particles. This will result in 
integrals being inversely proportional to the velocity of the particle. One solution is to 

20 normalize the area of each peak by the velocity of the corresponding particle to obtain 
velocity independent measurement of the fluorescence intensity. 

Several methods for measuring the velocity in microfluidic devices have 
been reported. In particle image velocimetry, video imaging is used to measure the 
velocities of particles in a channel by observing the displacement of the particles within a 

25 known time interval. See, for example, Singh et al.. Anal. Chem., 2001, 73, 1057-1061; 
Barker et al.. Anal. Chem., 2000, 72, 5925-5929; and Santiago et al.. Experiments in 
Fluids, 1998, 25, 316-319. This method is advantageous in obtaining the velocity spatial 
distribution, however it is not suitable for accurately measuring other analyte 
characteristic parameters, such as the fluorescent intensity. 

30 Shah convolution Fourier transform is another method to measure velocity 

in microfluidic devices. See, for example, Kwok et al.. Anal. Chem., 2001, 73, 1748- 
1753 and Crabtree et al., Anal. Chem., 1999, 71, 2130-2138. In this method, a mask with 
a periodic array of sUts spatially modulates the excitation beam. When an analyte is 
moving in the beam the spatial modulation is converted into a temporal modulation. The 



2 



distribution of velocities is found by Fourier transforming the temporal signal and 
identifying the peaks. However, such a practical use of this method has not been 
demonstrated. Moreover, the practical implementation of this method most likely will 
require fabrication of the mask on the chip which adds to the complexity of the device. 
5 Therefore, there is a need for apparatuses and methods for determining a 

velocity independent analyte characteristic parameter. 

SUMMARY OF THE INVENTION 
One aspect of the present invention provides an apparatus for determining 
a characteristic parameter of an analyte in a fluid medium independent of the flow 
10 velocity of the analyte, said apparatus comprising: 

(a) a device comprising a fluid flow channel; 

(b) a means for transporting a fluid medium from a first position to a 
second position of the fluid flow channel; 

(c) a plurality of detection zones located at different positions along 
1 5 said fluid flow channel and located in between the first and the second position of the 

fluid flow channel; 

(d) a detector for detecting the analyte flowing through the detection 

zone; and 

(e) a means for measuring a characteristic parameter of the analyte 
20 independent of the flow velocity of the analyte through the plurality of detection zones. 

Preferably, the device which comprises a fluid flow channel is a 
microfluidic device. While the fluid medium can be transported through the fluid flow 
channel by any conventional means, preferably means for transporting a fluid medium 
comprises a peristaltic pump or electroosmosis. More preferably, the microfluidic device 
25 comprises a peristaltic pump. Such microfluidic devices are particularly advantageous 
because the median fluid medium flow velocity can be controlled. 

In one particular embodiment, the detector comprises a laser, a laser beam 
guiding device and a fluorescence measuring device. Preferably, the laser beam guiding 
device is an acousto-optic modulator. In this manner, a velocity independent integrated 
30 fluorescence peak area of the analyte can be detemined. 

Another aspect of the present invention provides a method for determining 
a velocity independent characteristic parameter of an analyte, wherein the characteristic 
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parameter of the analyte is capable of being influenced by or dependent on the velocity of 
the analyte, said method comprising: 

(a) providing a means for transporting a fluid medium comprising the 
analyte from a first position to a second position of a fluid flow channel of a fluidic 

5 device; 

(b) measuring the characteristic parameter of the analyte within the 
fluid flow channel at a plurality of locations along the fluid flow channel in between the 
first and the second position; and 

(c) determining the velocity independent characteristic parameter of 
10 the analyte using the measured characteristic parameters of step (b) and normalizing the 

measurement by substantially eliminating the velocity component of the measurement. 

The velocity independent characteristic parameter of the analyte can be 

determined by: 

(i) comparing signals obtained from the plurality of locations along 
1 5 the fluid flow channel in step (b); 

(ii) determining a time difference by calculating the time it takes for a 
particular analyte to pass from a first detection position to a second detection position; 
and 

(iii) determining the velocity independent characteristic parameter 
20 using the time difference. 

Preferably, the signals from the first and the second detection zones are 
averaged and normalized using the time difference. By averaging the signals from the 
two detection zones, the amount of noise and false signals are significantly reduced by the 
methods of the present invention. 

25 Apparatuses and methods of the present invention can be used in a variety 

of assay and analytical apphcations. In one particular embodiment, methods of the 
present invention provide cell sorting. Yet in another embodiment, methods of the 
present invention provide determining number of nucleotides present in an 
ohgonucleotide. Such methods for determining the number of nucleotides in an 

30 oligonucleotide comprise: 

(A) attaching a fluorescent molecule to the oligonucleotide to produce 
a modified ohgonucleotide prior to measuring velocity independent characteristic 
parameter of the modified ohgonucleotide, wherein said characteristic parameter is 
integrated fluorescent peak area of said modified ohgonucleotide; and 
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(B) determining the number of nucleotides in the oligonucleotide by 
comparing the velocity independent integrated fluorescence peak area of the modified 
oligonucleotide with a velocity independent fluorescence peak area of a standard 
oUgonucleotide, wherein the velocity independent fluorescence peak area of the standard 
5 oligonucleotide has been calibrated to the number of nucleotides present the standard 
oUgonucleotide. 

BRIEF DESCRIPTION OF THE INVENTION 
Figure 1 A is a fluorescence signal illustration showing that at a similar 
flow velocity the peak intensity and the peak area is proportional to the length of the 
10 DNA being detected by flow cytometry. 

Figure IB is a fluorescence signal illustration showing peak areas of two 
similar length DNAs with different flow velocity. 

Figure 2 A shows a schematic illustration of an apparatus of the present 

invention. 

1 5 Figure 2B shows CCD camera image of fluorescein solution flowing 

through a T-chaimel of a microfluidic device illuminated by the scanning beam and two- 
dimensional laser beam scanning by acousto-optic modulator and the corresponding 
square wave (top, i.e., Chi) and fluorescence data (bottom, i.e., Ch2). 

Figiu-e 2C illustrates a laser beam location and the order after passing 
20 through two orthogonally located acousto-optic modulators. 

Figure 3 shows histogram of the areas obtained from the fluorescence of 
two kinds of fluorescence beads with and without velocity normalization. 

Figure 4 shows the velocity as a function of time of an experiment with the 
same two kinds of fluorescent beads as in Figure 3. 
25 Figure 5 shows histogram of the velocities of 0.2 [im fluorescent beads 

that were in the system (soUd line) and the corresponding theoretical one (dash line). 

Figure 6 shows histogram of A- DNA, where the center peak represents A, 
DNA (48 kbp) and the small peak on the right represents pairs of hybridized X (96 kbp). 

Figure 7 shows a flow chart of a computer program for recording square 
30 waves and fluorescent peak data. 

Figure 8 shows a flow chart of a computer program for determining a 
velocity independent flow cytometry. 
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DEFINITIONS 

"Anal3^e" refers to any material in a fluid medium which can be analyzed 
using a detector. Exemplary analytes include cells; oligonucleotides, such as DNA, RNA 
and PNAs; other organic compounds, such as pharmaceutically active compounds 
5 including antibiotics, antiviral compounds, anticancer compounds, etc.; beads; resins; 
polymers; and the like. 

The terms "character of an analyte," "character parameter of an analyte" 
and "an analyte's characteristic" are used interchangeably herein and refer to a physical or 
chemical characteristic parameter of the analyte. Exemplary analyte's characteristics 
10 include molecular size (e.g., weight or length), fluorescence, infrared or UV/VIS 
absorption, nuclear magnetic resonance (i.e., NMR) spectrum, cell type, and other 
characteristics which can be measured or detected by a suitable detector known to one 
skilled in the art. 

The terms "fluid" and "fluid medium" are used interchangeably herein and 
1 5 refer to a gas, or preferably liquid. 

The terms "velocity" and "flow velocity" when referring to an analyte are 
used interchangeably herein and refer to the flow velocity of the analyte within a fluid 
medium. 

DESCRIPTION OF THE SPECIFIC EMBODIMENTS 
20 The present invention provides an apparatus for determining a plurality of 

characteristic parameters of an analyte in a fluid medium. In particular, the apparatus of 
the present invention allows determination of an analyte characteristic parameter and the 
analyte flow velocity within a fluid flow channel of a device. Preferably, the analyte 
characteristic parameter is capable of being influenced and/or affected by the analyte flow 
25 velocity. In this manner, the apparatus of the present invention provides or can be used to 
determine a velocity independent analyte characteristic parameter. 

The analyte characteristic parameter can be any physical or chemical 
parameter that can be measured or determined using conventional detecting means. 
Suitable analyte characteristic parameters include UVA/IS absorption, fluorescence, 
30 nuclear magnetic resonance, infrared spectrum, and other physical or chemical parameters 
which can be measured as an integrated area under a peak or a curve. For the sake of 
brevity and clarity, the present invention will now be described in reference to measuring 
fluorescence of the analyte. Of course in order to measure fluorescence, the analyte must 
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be capable of fluorescing when exposed to an appropriate electromagnetic radiation, e.g., 
ultraviolet and/or visible light. Thus, in one aspect the analyte comprises a fluorescent 
moiety. 

The total fluorescence intensity integrated over the excitation duration is 
5 velocity dependent measurement. As expected, when fluorescence of anal3^es is 

measured in a narrow detection zone (i.e., area), analytes with faster flow velocity spend a 
relatively small amount of time within the detection zone resulting in a relatively smaller 
fluorescence intensity peak width compared to slower analytes. For example, as 
illustrated in Figure 1 A, at a similar velocity the peak intensity and the peak area is 

10 proportional to the length of the DNA being detected by flow cytometry. However, as 
shown in Figure IB, if two similar length of DNAs have different velocity, the faster 
moving DNA will have smaller peak area (in Figure IB, the dotted line shows that the 
peak heights of both particles are the same) and the slower moving DNA will have larger 
peak area. This velocity difference can lead to misleading or erroneous interpretation 

1 5 flow cytometry data. 

The apparatuses and methods of the present invention significantly reduce 
or ehminate the affect of molecule velocity through a detection zone on the peak area by 
ehminating or normalizing the velocity factor of the analyte characteristic parameter. In 
particular, apparatuses of the present invention comprise at least two different detection 

20 zones along the analyte's flow path (i.e., fluid flow channel) to determine the velocity of 
each analyte that flows through the detection zones. By placing two different detection 
zones at a predetermined distance (i.e., "d") from each other, one can measure the 
velocity of the analyte flowing through the detection zones by measuring the time 
difference (i.e., "t") at which the analyte passes through the first detection zone and the 

25 second detection zone. Since the velocity (i.e., "v") is distance divided by time, the flow 
velocity of analyte is calculated by the formula: v=d/t, where v, d and t are those defined 
above. The peak area is then multiplied by the velocity (or simply divided by time since 
d is constant) to normalize the peaks, i.e., to eliminate the velocity factor. In this manner, 
a more accurate determination of the analyte characteristic parameter can be made. 

30 Preferably, each analyte passes through the detection zone individually, 

i.e., separately. Therefore, the variables such as flow channel width, concentration of the 
analyte, and other variables which can affect the number of analyte passing through the 
detection zone is adjusted such that a single analyte flows through the detection zone at 
any given time. For measuring characteristic parameters of cells, the width of fluid flow 



7 



channel is preferably in the order of from about 1 |im to about 1000 ^im, more preferably 
from about 10 jxm to about 100 j^m, and still more preferably from about 5 jam to about 
50 jam. 

For example, for measuring characteristic parameters of a relatively large 
5 molecules having molecular weight of upto about a few hundred MDaltons, such as 

oHgonucleotides including DNAs, RNAs, PNAs and hybrids thereof (having upto about 
few hundred base pairs); peptides; polymers; and other organic compounds, the width of 
the fluid flow channel is preferably in the range of from about 1 |um to about 50 p,m, and 
more preferably from about 3 [im to about 10 |Lim. While the above fluid flow channel 
10 width are provided as being particularly suitable for certain analytes, it should be 

appreciated that other fluid flow channel widths are also within the scope of the present 
invention. 

Devices comprising such fluid flow channel dimensions are well known to 
one skilled in the art. For example, any micro fluidic devices currently known to one of 

1 5 ordinary skill in the art can be used in the present invention. However, preferred 

microfluidic devices are constructed of single or multilayer soft lithography (MLSL) as 
described by Unger et al. in Science, 2000, 288, 1 13-116, and fixrther detailed in 
commonly assigned U.S. Patent Application Serial No. 09/605,520, filed June 27, 2000, 
which are incorporated herein by reference in their entirety. Other preferred microfluidic 

20 devices are disclosed in a commonly assigned U.S. Patent Application entitled 

"Microfluidic Devices and Methods of Use," which is filed even date with the present 
application and is ftirther identified by attorney docket No. 020174-002510US, and is 
incorporated herein by reference in its entirety. Moreover, specific examples of 
microfluidic flow cytometry for sorting cells and DNA's are disclosed in commonly 

25 assigned U.S. Patent Application Serial No. 09/325,667 and the corresponding published 
PCT Patent Application No. US99/13050, and U.S. Patent Application Serial No. 
09/499,943, respectively, all of which are incorporated herein by reference in their 
entirety. 

While one can use multiple detectors and electromagnetic radiation 
30 sources (e.g., laser for laser induced fluorescence), it has been found by the present 

inventors that one or more, preferably two, acousto-optic modulators in conjunction with 
an aperture is particularly suitable for providing two different detection zones from a 
single laser source. An acousto-optic modulator is readily available from a variety of 



sources including at http://www.brimrose.coni/acousto_modulators.html, which also 
includes a general discussion on the theory behind acousto-optic modulators (Brimrose 
Corp., Baltimore, Maryland). Other devices which is capable of guiding the laser beam 
into two or more different positions can also be used instead of an acousto-optic 
5 modulator. Such devices are well known to one of ordinary skill in the art and include 
rotating mirrors, gratings and other electromagnetic wave diffracting devices. The 
aperture allows emission of only one particular diffracted beam to illuminate the detection 
zones and blocks other diffracted laser beam. 

When using an acousto-optic modulator, preferably the first order beam is 

10 used. Thus, when two acousto-optic modulators are used, one to control direction of the 
laser beam in the x-axis and the other to control direction of the laser beam in the y-axis, 
the resulting laser beam is about 20% intensity of the original laser beam. This is because 
a typical first order laser beam intensity exiting an acousto-optic modulator is about 50% 
of the laser beam entering it. Therefore, by using two acousto-optic modulator, the 

1 5 resulting laser beam is about 20% to about 30% intensity of the original laser beam. 

As stated above, an aperture is typically used in conjunction with acousto- 
optic modulators. The aperture is used to select the first order laser beam firom each of 
the acousto-optic modulators. For example, as shown in Figure 2C, there are variety 
orders of diffracted laser beam resulting form two acousto-optic modulators. The 

20 aperture allows laser beam from only the first order of both acousto-optic modulator (i.e., 
50) to be focused onto the detection zones. 

A typical set-up for using a laser to detect fluorescence of the analyte in 
two detection zones is schematically illustrated in Figure 2A. The laser beam enters two 
orthogonally located acousto-optic modulators 100 A and lOOB. The aperture 104 then 

25 allows only the first order laser beam from both of the acousto-optic modulators to pass 
through. The resulting laser beam then passes through two lenses 108A and 108B. 
However, the use of two lenses is optional and can be omitted. Typically, the lenses are 
used to adjust the distance between the two scan lines (i.e., detection zones) as well as 
their length, for example, to be in the order of about 10 }xm. 

30 The laser beam is then deflects off the dichroic beam-splitter 112 and is 

focused on to the detection zones through an objective 116. The fluorescence wavelength 
then passes through a filter 120 to a channel photo-multipher 124, which is operatively 
interconnected to a data acquisition device 128. The optional oscilloscope 132 can be 
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used to observe the signal prior to data acquisition and/or to observe the signal in real- 
time. 

The result of one particular embodiment of the present invention is shown 
in Figure 2B. The left portion of Figure 2B shows a T-shaped flow channel which is 
5 illuminated with fluorescence at two different detection zone of the flow channel. The 
right portion of Figure 2B shows signal peaks corresponding to the square wave which 
controls the x-axis position, i.e., the detection zone, of the laser beam and the 
fluorescence peaks obtained from the two detection zones. In this embodiment, the 
distance 10 between the two detection zones 14L and 14R is about 10 jxm and the width 

10 of the fluid flow channel is about 5 ^im. The laser beam enters two orthogonally 

positioned acousto-optic modulators lOOA and lOOB and the resulting first order beam is 
emitted through an aperture 104. In order to scan the entire cross section of the fluid flow 
channel and to allow scan of two different detection zones 14L and 14R, two acousto 
optic modulators are used. One to control the x-axis and the other to control the y-axis of 

15 the laser beam. In Figure 2B, the beam has y-axis frequency of 150 kHz, i.e., the beam 
travels from above the "top" of the flow channel 20 to below the "bottom" of the flow 
channel 24 at a rate of 1 50,000 times per second. In actuality the beam scans twice the 
distance relative to the width of the fluid flow channel; however, because there is no 
fluorescence outside the fluid flow channel, no beam is visible outside the fluid flow 

20 channel in Figure 2B. Furthermore, the laser beam switches from the detection zone 14L 
to 14R and vice a versa at a rate of 5 kHz. In addition, the laser beam has a sampling rate 
of 40 kHz, i.e., each x-position is sampled about 4 times or 8 times total (40/5). 
Frequency of x-axis switching can be seen in the top graph of the right portion of Figure 
2B. In this graph, when the peak is at the top, it represents detection (or scanning) in the 

25 14R region, and when the peak is at the bottom (i.e., 0) it represents detection (or 

scanning) in the 14L region. As can be seen, the laser beam moves from one position to 
another (in the x-axis) to allow scanning of two different positions. This allows the same 
analyte to be detected at two different times at two different regions as shown in the lower 
graph of the right portion of Figure 2B. By determining the time difference between such 

30 detection and knowing the distance 10 (Figure 2 A), one can calculate the velocity of the 
material traveling through the fluid flow channel. As stated above, it is preferred that 
statistically each analyte enters the detection zone separately. 
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As stated above, the left portion of Figure 2B illustrates scanning two 
detection zones of a T-shaped fluid flow channel by the scanning beam. A fluid medium 
comprising fluorescein was introduced to the fluid flow channel and was excited by the 
scanning beam in two scan lines as shown. The resulting fluorescence was imaged with a 
CCD camera as shown m Figure 2B. The two scan lines are clearly visible. It is 
important to note that the two scan lines can be seen together only due to the limited time 
resolution of the CCD compared to the scanning frequency. In actuality, the two bright 
lines of Figvtre 2B are actually fluorescing at different times. The lines scanned by the 
illuminating beams are about two times longer than the bright lines shown in Figure 2B to 
ensure uniform excitation of the region of interest. However, the full line can not be seen 
because there is no fluorescein outside of the fluid flow channels, and therefore no 
fluorescence occurs there. The width of the fluid flow channel in Figure 2B is 5 |am and 
its depth is 3.9 |am. 

As shown in Figure 2 A, in one embodiment the laser beam is transmitted 
through two mutually orthogonal acousto-optic (AO) modulators, represented as items 
lOOA and lOOB, to produce two line scans. Each AO was driven by voltage controlled 
oscillator (VCO) connected to a function generator (not shown), such as Stanford 
Research Systems, Model DS335 (Sunnyvale, CA). As implied above, one of the AO is 
used to move the beam in the x direction, i.e., parallel to the direction of the fluid flow. 
Typically, this VCO was driven with a square wave (e.g., 5KHz, 7Vp-p) thereby allowing 
one to scan two different detection zones along the fluid flow channel. Because the laser 
beam scans only a small fraction of the width of the fluid flow channel, a second AO is 
required to scan the entire width of the fluid flow channel. This second AO is preferably 
oriented to enable the laser beam to move in the Y direction, thus allowing the beam to 
cross the entire width of the fluid flow channel. Typically, the signal for controlhng this 
second AO is a sinusoidal wave (e.g., 120KHz, 9Vp-p) at a frequency higher than the 
chopping signal to ensure several crossing of the channel for each position on the x 
direction. At the exit of the second AO the original laser beam is divided into a grid of 
beams as illustrated in Figure 2C. Most of these beams are blocked using an aperture 104 
leaving only the beam that is the combination of the first orders of both AOs. This beam 
is then used to scan along two parallel lines (left and right line scans) in different 
positions along the fluid flow. Such scanning can be conveniently accomplished by a 
VCO. As stated above, the data acquisition can be conveniently accomphshed by a 
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channel photo-multiplier 124 operatively connected to a data acquisition device, such as a 
computer 128 or other signal recording device, including a plotter. Typically, the two 
signals, the square wave used to control AOs and the fluorescence signal, are recorded 
separately, i.e., as Chi and Ch2 in Figure 2B. By correlating signals of Cli2 to the square 
5 wave of Chi allows one to determine which peaks in Ch2 are from the left (14L) 
detection zone and which peaks are from the right (14R) detection zone. 

Figure 3 is a histogram of the peak areas for 14L and 14R regions, as well 
as normalized, i.e., velocity independent, peak areas for an experiment conducted with 
two kinds of beads which differ in their fluorescent intensity. It is expected that the ideal 

10 histogram will have two peaks, one for each kind of beads. As the legends in Figure 3 
shows, the graph represents none normalized fluorescence peak areas from the left (14L) 
and right (14R) detection zones and a normalized, i.e., velocity independent, peaks 
(labeled as "both"). It has been found that in some cases averaging the normalized peak 
areas can reduce the coefficient of variance relative to using only a single detection zone. 

15 It has been found by the present inventors that because apparatuses of the 

present invention provide a pliu-ality of detection zones, the resulting peak area (i.e., 
characteristic parameter) that has been normalized have a significantly lower coefficient 
of variance compared to conventional apparatuses having only one detection zone. 
Typically, apparatuses and methods of the present invention improves the coefficient of 

20 variance by a factor of at least about 2 relative to a similarly equipped apparatus having 
only one detection zone, and preferably by a factor of at least about 3. 

Gel-electrophoresis (i.e., electrophoresis) and other similar methods have 
limited resolution capacity for medium to large DNA molecules, and therefore are 
inapplicable in many cases. In contrast, apparatuses and methods of the present invention 

25 are not limited by the size of material (e.g., DNA). Moreover, if the distance between 
two detection zones are large or the velocity of the material is slow, one can use these 
variations to study a variety of analyte characteristic parameters. For example, one can 
detect changes in cells as it passes through from one detector to another. One can also 
analyze chromosome distribution in cells (e.g., karyotyping). Methods of the present 

30 invention are also useful in epidemiology and other diagnostic and assay procediu-es. 

In one particular aspect, apparatuses of the present invention can be used 
to determine the size (e.g., number of base pairs) of oligonucleotides such as DNAs, 
RNAs, PNAs and hybrids thereof. For example, a calibration chart can be prepared by 
measuring fluorescence of oligonucleotides of a variety of lengths. To be useftil, each 
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oligonucleotides are coupled to a compound that is capable of fluorescing. The 
fluorescence compound is then coupled to the ohgonucleotide at a regular interval, e.g., 
every four or five base pairs. The oligonucleotide of unknown length is then coupled to 
the same fluorescence compound at the same nucleotide interval. By measuring the 
5 fluorescence peak area and comparing the result with the calibration chart, e.g., by a 
computer, one can easily determine the length of the unknown ohgonucleotide. Thus, 
apparatuses and methods of the present invention can be used as an alternative to 
electrophoresis to determine the size of ohgonucleotides. However, imlike 
electrophoresis, apparatuses and methods of the present invention are generally not 

1 0 limited by the size of the ohgonucleotide. 

Apparatuses and methods of the present invention allow determination of 
analyte flow velocity and other characteristic parameter(s) of the analyte. For example, 
such apparatuses and method can be used to perform flow cytometry and determine the 
anal3^e flow velocity; thereby, enabhng correlation of the velocity with each specific 

1 5 analyte. In addition, usefiilness of apparatuses of the present invention is not limited by 
the size of the analyte. In contrast, measuring the fluorescence fi:om the peak heights 
using a conventional apparatus typically requires the entire particle (i.e., analyte) to be 
uniformly illuminated by the excitation beam, thereby Umiting the size of the particle 
which can be analyzed. 

20 Data Acquisition and Analvsis 

In one particular embodiment of the present invention, the data acquisition 
involves scanning two different detection zones along the path of the fluid flow channel 
by AOs which are operatively interconnected to VCOs. Each detection zone is scanned at 
least the entire width of the fluid flow channel. Typically, each acquisition involves ten 

25 or more scans per detection zone. That is, the laser beam scans "up and down" the width 
of the detection zone (i.e., y-axis of the fluid flow channel) ten or more times before the 
X-axis position (i.e., detection zone) of the beam is switched. Data acquisition typically 
involves recording the fluorescence signal in each of the two detection zones and 
separately recording the square wave that is used to direct the laser beam from to different 

30 detection zones. If the peak and the valley of the square wave represent relative values 1 
and 0, respectively, typically 0.05 and 0.95 firactions of square wave signals are used as 
the locations for the first and the second detection zones. This eliminates any noise factor 
that may be present in the signal, e.g., due to signal spikes. 
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As stated above, data acquisition involves recording the square wave 
frequency and the fluorescence intensity of the analyte. By correlating the square wave to 
the fluorescence peaks, one can determine whether a particular fluorescence peak is from 
the first or the second detection zone. For example, if the laser beam is directed at the 
first detection zone while the square wave is near its peak value, then any fluorescence 
peaks occurring at the same time as the peak of the square wave is due to fluorescence at 
the first detection zone. Conversely, any fluorescence peaks occurring while the square 
wave is near its bottom is due to fluorescence at the second detection zone. In this 
manner, the fluorescence signal can be separated into its two origin signals (one for each 
detection zone) even in cases where a particle is long enough to be in both detection 
zones at the same time. This type of data analysis algorithm allows detection of 
fluorescence from both detection zones using a single detector. 

After each peak from the first detection zone has been correlated to a 
corresponding peak from the second detection zone, the peaks are then normahzed by 
dividing the integrated peak area with the analyte flow velocity. The analyte flow 
velocity can be determined by the time difference between the positions of the two 
corresponding peaks. Since the distance between the first and the second detection zone 
is constant, the integrated peak area can be simply divided by the time it take for the 
analyte to travel from the first detection zone to the second detection zone. Because each 
integrated peak area from the first detection zone is correlated to the corresponding 
integrated peak area from the second detection zone, any integrated peak area from the 
first detection zone that does not have a corresponding integrated peak area in the second 
detection zone or vice a versa is most likely due to noise and is discarded. Thus, 
apparatuses and methods of the present invention provide more accurate analyte 
characteristic parameter determination than conventional single detection methods. 

The distance between the two detection zones should be sufficiently large 
enough to allow accurate velocity determination, i.e., time it takes for the analyte to travel 
from the first detection zone to the second detection zone. Generally, the accuracy of the 
velocity determination is somewhat dependent on separation (in time or more importantly 
in sampling points) between the two corresponding fluorescence peaks. Therefore, larger 
separation of the fluorescence peaks, i.e., longer it take for an analyte to travel from the 
first detection zone to the second detection zone, results in more accurate velocity 
determination. 
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In one particularly preferred embodiment, statistically only one analyte 
passes through both detection zones before the second particle enters the first detection 
zone. In this manner, the integrated peak areas from the first detection zone can be 
readily correlated to the integrated peak areas from the second detection zone. This can 
5 be accomplished by adjusting one or both of the distance between two detection zones 
and the concentration of the analyte in the fluid medium. Of course these two parameters 
are also dependent on other variables such as the mean flow velocity of the fluid medium 
within the fluid flow channel and the square wave frequency (i.e., frequency at which the 
laser beam switches to and from the first and the second detection zones), 

1 0 For example, for a mean fluid medium flow velocity of about 2 lam/sec 

and square wave frequency of 5000 sec'^ the distance between the two detection zone is 
about 10 j^m. And the concentration of the analyte is typically from about 2 nM or less, 
preferably from about 0.2 nM to about 2 nM. 

Theoretically, activating data acquisition prior to having any of the analyte 

15 pass through the first detection zone allows one to easily correlate each fluorescence peak 
to a particular particle, becasue the order of each peak on the first and the second 
detection zones will be equal to the order of particles passing through each detection 
zone. However, in many instances such signal acquisition is not feasible or practical. 
Moreover, possible noises in each of the detection zones will prevent correlating the peak 

20 from one detection zone to another detection zone. Therefore, one needs to correlate the 
integrated peak area (e.g., fluorescence peaks) from the first detection zone to the 
corresponding integrated peak area of the second detection zone. 

One method of achieving this correlation is to correlate the integrated peak 
area from the first detection zone to a corresponding integrated peak area in the second 

25 detection zone by limiting the detection to region of velocities about ± 4 times the median 
flow velocity of the fluid medium. For example, integrated peaks from the first detection 
zone is compared to integrated peaks from the second detection zones in the range of 
the median flow velocity to 4x the median flow velocity of the fluid medium. By 
comparing other integrated peak areas of similar intervals, one can verify whether the 

30 initial correlation is accurate. Such a method for detecting and correlating the integrated 
peak areas from the first detection zone to the integrated peak areas of the second 
detection zone is schematically outlined in the flow sheets in Figures 7 and 8. Briefly, the 
two integrated peak areas (e.g., Ai(i) and Ar(i)) are compared to see whether they are 
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within the minimum and the maximum time difference as specified. The minimum and 
maximum time differences can be adjusted depending on the median flow velocity of the 
fluid medium. Once a pair of matching integrated peaks from the first and the second 
detection zones is found, the time difference (i.e., TimeDiff or At) is determined. The 
5 time difference is inversely proportional to the analyte flow velocity and is the difference 
in time when the particle crosses (i.e., detected by) the first and second detection zones. 
The velocity independent integrated peak area is then calculated by averaging the 
integrated peak areas of the particle from the first and the second detection zones and 
dividing the average integrated peak area by time difference. This calculation can be 
1 0 represented by the following formula: 

Avi = [(A] + Ar)/2]/At 

where Avi is velocity indpendent area, Ai is integrated peak area from the first detection 
zone, Ar is integrated peak area from the second detection zone and At is time difference. 
In Figure 8, the time difference is determined by dividing the difference in point number 

15 of A] and Ar with the sampling rate (i.e., the frequency of switching the laser beam to and 
from the first and the second detection zones). 

Since the velocity independent peak area is determined by averaging the 
two integrated peak areas from the first and the second detection zones, the apparatuses 
and methods of the present invention reduces the noise significantly. Moreover, any 

20 integrated peak from one detection zone that does not have a corresponding integrated 
peak on the second detection zone or vice versa is eliminated resulting in further 
reduction of noise. In conventional single detection method, such peaks are often 
mistaken as an actual signal, thereby resulting in misinterpretation of the data. By 
comparing the signals from one detection zone to the other, methods of the present 

25 invention further reduces noise or false signals. 

It should be appreciated that while the present invention has been 
described with respect to using a single laser beam with two acousto-optic modulators to 
guide the laser beam into two scanning detection zones, other arrangements are also 
possible and are within the scope of the present invention. For example, apparatus of the 

30 present invention can be fabricated to utilize two laser beams directed at different 

detection zone and having a separate acousto-optic modulator for scanning the entire 
width of the fluid flow channel. In addition, other electromagnetic wave diffracting 
devices, such as rotating mirrors, gratings and other electromagnetic wave diffracting 
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devices known to one skilled in the art, can be used instead of an acousto-optic 
modulator. Other variations include using two detectors with either one or two laser 
beams. If one laser beam is used, it can be spUt or scanned (i.e., guided) at two different 
position as described herein, for example, with AO, mirror or other suitable means. If 
two lasers are used, no scanning is necessary as each laser can be used to illuminate 
different detecting area. These setup variations offer easier data acquisition and analysis; 
however, in some cases may result in more comphcated aUgnment of each components. 

Moreover, the control over the line scans (i.e., detection zones), which is 
done by controlling the AOs, is done electronically and can be changed depending on a 
particular application. For example, the two line scans can be easily rotate by 90 degrees 
to enable the same measurement in a perpendicular channel. This general measurement 
method can be used to measure the velocity of particles in a variety of microfluidic 
devices without the need for any changes in the design of microfluidic devices. 

Microfluidic Device 

Apparatuses of the present invention also comprise a fluidic device, 
preferably a microfluidic device, which comprises a fluid flow channel. The fluid flow 
channel allows flows of the fluid medium from one location to another location within the 
device. Preferably, the cross-section of the fluid flow channel should be small enough 
such that only a very small area is scanned by a laser. In general, any microfluidic device 
made of a material that is transparent to the laser beam is suitable. However, a 
particularly preferred microfluidic devices are those disclosed by Unger et al. in Science, 
2000, 288, 1 13-1 16, and U.S. Patent Application Serial No. 09/605,520, which were 
incorporated by reference above. 

In particular, microfluidic devices which comprise a peristaltic pump is 
particular useful in the present invention as they allow one to control the median fluid 
medium flow velocity. However, it should be appreciated that the fluid medium can be 
made to flow through the fluid chaimel by any of the conventional means, such as 
pressure gradient, electroosmotic flow, and the like. 

Utilitv 

Apparatuses and methods of the present invention have a wide variety of 
appUcation such as flow cytometry, ohgonucleotide sorting, oligonucleotide analysis, 
detecting changes in cells, analyzing chromosome distribution in cells (e.g., karyotyping), 
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studying epidemiology, and other diagnostic and assay procedures. Using microfluidic 
devices in cell and DNA sorting are general described by Fu et al. in "A Microfabricated 
Fluorescence-activated Cell Sorter " Nature Biotech., 1999, i 7, 1 109-1 1 1 1 ; and Chou et 
al. in "A Microfabricated Device for Sizing and Sorting DNA Molecules," Proc. Natl 
5 Acad. Sci. USA, 1999, 96, 11-13, respectively, which are incorporated herein by reference 
in their entirety. Methods for cell sorting using microfluidic devices are fiirther detailed 
in PCT Publication No. WO 99/61 888, which is also incorporated herein by reference in 
its entirety. 

In one particular embodiment, methods of the present invention provides 
1 0 determining number of nucleotides that is present in an oligonucleotide. The method 
generally involves modifying the oligonucleotide by attaching a fluorescent molecule to 
the oligonucleotide and measuring a velocity independent integrated fluorescent peak 
area. This velocity independent integrated fluorescent peak area is then compared with a 
velocity independent integrated fluorescent peak area of a standard oligonucleotide 
1 5 having a known number of nucleotides. Preferably, velocity independent integrated 

fluorescent peak areas are determined for a number of standard oligonucleotides having a 
different number of nucleotides. In this manner, a calibration table can be produced and 
used to determine the number of nucleotides present in an oligonucleotide. 

The modified oligonucleotide and the standard oligonucleotides are 
20 attached with a known amount of fluorescent molecules per given number of nucleotides. 
In this manner, a given velocity independent integrated fluorescent peak area can be 
calibrated to represent a certain nimiber of nucleotides that is present in the 
ohgonucleotide. 

Additional objects, advantages, and novel features of this invention will 
25 become apparent to those skilled in <he art upon examination of the following examples 
thereof, which are not intended to be limiting. 

EXPERIMENTAL 

Apparatus 

The fluorescence excitation was performed with a 5 mW, 488 nm laser 
30 beam fi-om an air cooled argon ion laser (Uniphase, San Jose, CA). The laser beam was 
transmitted through two mutually orthogonal acousto-optic (AO) modulators to produce 
two Hne scans (see Figure 2A). Each AO was driven by voltage controlled oscillator 
(VCO) connected to a function generator. The first AO was oriented so as to move the 



beam in the x direction, parallel to the direction of the flow. Its VCO was driven with a 
square wave (5 KHz, 7 Vp.p) causing the beam to move between two points along the 
channel. The second AO was oriented to enable moving of the beam in the Y direction, 
thus allowing the beam to cross the channel. The signal that was used to control this AO 
5 was a sinusoidal wave (120 KHz, 9 Vp.p) at a frequency higher than the chopping signal 
to ensure several crossing of the channel for each position on the x direction. At the exit 
of the second AO the original laser beam was divided into a grid of beams (see Figure 
2C). Most of these beams were blocked leaving out only the beam that was the 
combination of the first orders of both AOs. This beam was thus scanned along two 

10 parallel lines (left and right line scans) in different positions along the flow. The laser 

beam was focused through a 100 X 1.3NA oil immersion objective (Olympus, New Hyde 
Park, NJ) which also was used to collect the emitted fluorescence. The full width half 
maximum of the beam at the channel was about 1 \xm. Auxihary lenses were used to 
adjust the distance between the two scan lines as well as their length to be in the order of 

15 10 |j,m. The uniformity of the excitation across the channel was not dependent on the 
beam size but rather on the uniformity of the illumination while the beam was scanned. 
This was evaluated by imaging a thin layer of fluorescein in solution illuminated by the 
scanning beam with a CCD camera (see Figure 2B, left portion). The image was then 
digitized and evaluated for uniformity. A dichroic filter was used to introduce the laser 

20 hght into the optical train (Chroma 500 DCLP, Chroma Technology, Brattleboro, VT). 
Dielectric filter was used to reduce background and scattered light from the emitted 
fluorescence (Chroma D535/50M). The fluorescence was imaged onto a channel 
photomultiplier detector (EG&G, Gaithersburg, MD). The detector output as well as the 
chopping square wave were digitized at a rate of 40 KHz by a National Instrument 

25 (Austin, TX) Lab PC 1 200 board on a personal computer running LAB VIEW. 
Microfluidic device 

A 3.9 jxm high, 5 \xm wide rectangular fluid flow channel was fabricated 
from a silicon elastomer (General Electric RTV 615) by using a replica technique 
generally disclosed by Chou et al., in Proc. Natl. Acad. Sci. USA, 1999, 96, 11-13 and 

30 U.S. Patent Apphcation Serial No. 09/605,520, filed June 27, 2000, which were 

previously incorporated by reference in their entirety. Master molds were made from 
silicon wafers by using standard micromachining techniques. 
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Data Analysis 



Both line scans were well inside the detecting area, thus the fluorescence 



signals from both Une scans were recorded as one signal. However, since the excitation 
beam was scanning, the obtained fluorescence light originated from only one line scan. 
5 Thus, by simultaneously recording the data signal and the chopping square signal, it was 
possible to separate the data signal to two separate subsignals, one from each line scan 
(see flow chart in Figure 7). Each of these signals (5 KHz) was individually analyzed off 
line with a computer program to obtain the position in time and area of each peak in the 
signals. 

10 As shown in the flow chart in Figure 8, the positions of the peaks from 

both subsignals were used to match pair of peaks, one from each line scan, that were 
assumed to arise from the same fluorescent molecule. Then for each pair the time 
difference between the positions of its two peaks were calculated. The time difference is 
inversely proportional to the velocity of the molecule. The average area of each pair of 

1 5 peaks was normalized by dividing it by the corresponding time difference. Thus, for each 
pair of peak, which arises from one particle, the normalized total fluorescence and the 
velocity were obtained. 

Example 1 



0.2 \im., yellow green fluorescent. Molecular Probes) were introduce into the microfluidic 
device and their velocity were measured. The flow of fluorescent beads was achieved by 
capillary flow in the microfabricated channels sealed by coverslips. A histogram of the 
25 bead flow velocities is seen in Figure 5 for a later time in the experiment in which the 

transient stage is over and the velocity distribution become roughly constant. Assuming 
imiform distribution of beads in the channel, the theoretical probability, Ppr, of a particle 
with velocity vo to pass through the cross section in a given time is given by the equation: 



30 where Py is the probability to have a velocity v in the cross section of the chaimel. In 

circular tubes is constant for v<Vmax for Poiseuille flow. The rectangular cross section 
of the fluid flow channels of the microfluidic device had an aspect (i.e., width to length) 



This example illustrates the method of measuring flow velocity of 



20 particles. 



A solution of fluorescent beads (FluoSpheres Biotin Labeled microspheres 




Vo-Pv(v-Vo) 
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ratio of 1 .3; therefore, only minor changes to Py is introduced compare to a circular cross 
section fluid flow channel. Therefore, as a first order approximation the theoretical Ppr is 
proportional to v for v<Vmax- This theoretical result is substantially in agreement with the 
experimental results. See Figure 5. The high correlation between the two graphs shows 
5 that velocity can be accurately determined using an apparatus and a method of the present 
invention. 
Example 2 

In the next experimental stage a solution of a mix of two kinds of beads 
(Component A and B, LinearFlow green flow cytometry intensity calibration kit, 2.5 |am, 

10 Molecular Probes) which were similar m every other aspect besides their fluorescent 

intensity was introduce to the system. The fluorescent intensity of the beads was collected 
by the channel photomultipUer for a period of 400 sec. For each passing particle, the 
normahzed area and velocity were calculated. 

The velocity as a function of time was plotted. See Figure 4. The flow in 

15 the system is capillary flow and is driven by the pressure difference between the ends of 
the channel. During the time of the experiment, the pressure at the begirming of the 
channel, originating from the fluid at the starting basin, decreases while the pressure at 
the output increases. This press\u-e change results in a decay in the velocity over time, as 
shown in Figure 4 (t > 20sec). The abrupt change in the velocity, which occurs at about 

20 20 sec, corresponds to the arrival of the fluid front to the output end of the fluid flow 
channel. 

A histogram of the peak areas and the normalized peak areas was plotted. 
As shown in Figure 3, without normahzation it is difficult to distinguish between the two 
kinds of beads. The coefficient of variance (CV) of the peaks were improved by the 
25 normahzation from about 24% and 18% to 9.75% and 5.18% for the two kinds of beads, 
respectively. 

An extra peak in the histogram is seen only in the graphs from the single 
line scan at the lower area region, which is believed to be associated with noise. Many of 
the noise peaks that were found in each line scans signal can not be paired with 
30 correlating peaks from the other channel. Thus, the additional low peak in the histogram 
is missing for the normalized areas one, which shows that the normalization method 
reduces noise significantly. 



21 



Example 3 

Lambda phage DNA (GIBCO) was diluted in buffer (Tris EDTA, pH 6.8 
with 5 mM NaCl) and stained with the intercalating dye YOYO-1 (Molecular Probes) at a 
stochiometry of one dye molecule per 4 bp. 
5 Then a solution of stained X DNA was introduced to the microfluidic 

device for a period of 5 minutes. A histogram of the normalized areas of the molecules is 
plotted and is shown in Figure 6. The center peak in the histogram corresponds with the 
X DNA (48 kbp) and the small peak in the right corresponds with hybridized pairs of 
X DNA (96 kbp), i.e., X^. The CV for the X and X^ peaks were 8.34% and 7.61%, 
10 respectively, and the ratio between the centers of the two peaks was 2.06. 

The foregoing discussion of the invention has been presented for purposes 
of illustration and description. The foregoing is not intended to Umit the invention to the 
form or forms disclosed herein. Although the description of the invention has included 
description of one or more embodiments and certain variations and modifications, other 

1 5 variations and modifications are within the scope of the invention, e.g., as may be within 
the skill and knowledge of those in the art, after vmderstanding the present disclosure. It 
is intended to obtain rights which include alternative embodiments to the extent 
permitted, including alternate, interchangeable and/or equivalent structures, functions, 
ranges or steps to those claimed, whether or not such alternate, interchangeable and/or 

20 equivalent structures, functions, ranges or steps are disclosed herein, and without 
intending to publicly dedicate any patentable subject matter. 
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